Loop models and their critical points 
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Loop models have been widely studied in physics and mathematics, in problems ranging from 
polymers to topological quantum computation to Schramm-Loewner evolution. I present new loop 
models which have critical points described by conformal field theories. Examples include both fully 
packed and dilute loop models with critical points described by the superconformal minimal models 
and the SU "(2)2 WZW models. The dilute loop models are generalized to include SU(2) k models 
as well. 



I. INTRODUCTION 

The statistical mechanics of loop models have been 
studied extensively for decades. These models describe 
physical systems, and also provide fundamental problems 
in mathematics. Major progress has been made recently 
in understanding their properties. This is especially true 
when the loops are embedded in two dimensions, where 
the extensive variety of theoretical tools available have 
made many computations not only exact, but rigorous 
as well. 

In a loop model, the degrees of freedom are one- 
dimensional. These "loops" may branch and touch, but 
are not allowed to have ends. The partition function is 
of the form 
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no 
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This sum is made precise by defining the loops to live on 
the links of some lattice. Then £ labels a single loop con- 
figuration where the loops have total length L(C), i.e. L 
links of the lattice are covered by a loop. The parameter 
t is therefore a weight per unit length, which in this paper 
will usually be tuned to a special value which makes the 
behavior critical. Different loop models are distinguished 
by the choice of the weight w(C). Here it is required that 
w(C) depends only on topological properties of a given 
configuration C. Topological weights can depend on the 
number of loops, how many times loops touch, and more 
intricate properties like the number of ways of coloring a 
given graph. 

A simple and famous example of a loop model goes 
under a variety of names: "ring polymers", "closed self- 
avoiding random walks", or the "O(n) loop model" 
Here the loops are defined so that they do not branch and 
do not touch each other, i.e. if two or more loops touch 
the same point, a prescription for resolving them must 
be given. A typical configuration is displayed in figure [TJ 
the lattice is rotated by 45 degrees for later convenience. 
A simple topologically-invariant quantity is the number 
of loops Af in a given loop configuration C. The 0(n) 
loop model is then defined by taking 



w(C) = n M 



(2) 



FIG. 1: A typical configuration in the 0(n) loop model on 
the square lattice; the sites drawn are on the dual lattice. 



for some parameter n. A single loop/walk/polymer can 
be obtained by requiring n — > 0. An interesting varia- 
tion on this model comes from requiring that every site 
of the lattice be visited by at least one loop. Such loops 
are called fully packed, and their critical points are very 
different from those of the 0(n) model without this re- 
quirement [H, H, [H . 

Many well-known models of statistical mechanics can 
be expressed in terms of loops. A simple example is the 
Ising model in two dimensions. A loop configuration is 
found from each spin configuration simply by writing the 
domain walls separating regions of up and down spins. 
A given loop configuration corresponds to two spin con- 
figurations, since flipping all the spins does not change 
the domain walls. These Ising loops/domain walls can 
touch, but do not have ends, and cannot branch: there 
must be even number of links with loops at any lattice 
site. The weight per unit length t = e~ 2K for the usual 
Ising nearest- neighbor coupling K = J/(ksT), while the 
topological weight w(C) — 1 for all allowed C. 

It is essential to note that statistical-mechanical mod- 
els defined by |T]) do not necessarily have local Boltzmann 
weights. Since a single loop can be arbitrarily large, 
topological properties like the number of loops are de- 
fined non-locally. In some cases, there is a model with 
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local Boltzmann weights with the same partition func- 
tion. Generically, these Boltzmann weights are complex, 
and so the theory is non-unitary. In some special cases 
like the 0(n) model with n = 2cos(7r/(/c + 2)) with k a 
positive integer, a local model with positive Boltzmann 
weights does exist. 

Interesting questions in both physics and mathematics 
arise in the continuum limit. In this limit one can apply 
a variety of powerful field-theory methods to understand 
the phase diagram and compute properties at and near 
critical points. One major breakthrough in the study of 
two-dimensional loop models came with the invention of 
Coulomb-gas methods (see [l| for a review) . These meth- 
ods allowed non-trivial critical points to be found, and 
many critical exponents to be computed exactly, under 
a well-motivated set of assumptions about the behavior 
of the loops in the continuum limit. Further progress 
in understanding these critical points came with the in- 
troduction of conformal field theory 4]. When a critical 
point in a loop model is in the same universality class as 
a known conformal field theory, then not only can criti- 
cal exponents be found, but correlators computed exactly 

A remarkable recent development in the study of loop 
models came with the introduction of Schramm-Loewner 
evolution (SLE) [6]. SLE provides a novel way of for- 
mulating conformal field theory geometrically. Among 
many other things this gives a direct way of deriving 
which conformal field theory describes certain loop mod- 
els; sometimes the correspondence can even be proved 
rigorously. 

Another reason for recent interest in loop models comes 
from the search for quantum theories with quasiparticles 
obeying non-abelian statistics. Aside from their intrinsic 
interest, such quasiparticles could make up the "qubits" 
of a topological quantum computer 0, Q. One of the 
few types of models exhibiting such quasiparticles are 
quantum theories whose Hilbert space has basis elements 
consisting of loops in two dimensions. The ground-state 
wave function is a sum over all such loops, so that ground- 
state correlators in the quantum theory are identical to 
those in a two-dimensional classical loop model like the 
ones studied here. 

Despite all the recent activity and progress, not many 
different types of loop models have been studied in depth. 
One reason is that many loop models do not have critical 
points. For example, the 0(n) loop model has a critical 
point only when n < 2. There are many interesting and 
simple conformal field theories which have no known de- 
scription in terms of loops. Moreover, for all the successes 
of the SLE approach, it applies only to a very small set 
of conformal field theories. 

The purpose of this paper is to develop a general pro- 
cedure for finding a two-dimensional loop model with a 
critical point and a corresponding conformal field theory 
description. This procedure yields a number of new mod- 
els whose geometric properties should be understandable 
by using conformal field theory, and will hopefully prove 



useful for SLE and for topological quantum computation. 
For example, the loop models described here may provide 
a way of geometrically interpreting the SLE processes de- 
scribed in 0]. 

The strategy is to start with classical lattice models 
whose degrees of freedom are not loops, but rather are 
heights, integer- valued variables on the sites of the square 
lattice. The Boltzmann weights are local and positive, 
so any field theories obtained from taking the continuum 
limit of these models are unitary. The height models uti- 
lized are integrable, which makes it possible to find a va- 
riety of exact results. Choosing these models judiciously 
makes it possible to 

1. Locate critical points within the lattice model, and 
find the conformal field theory describing the con- 
tinuum behavior of each critical point. 

2. Re-express the partition function in terms of loop 
degrees of freedom, instead of heights, just as dis- 
cussed for the Ising model above. 

Putting these two together gives a loop model whose crit- 
ical point is described by a known conformal field the- 
ory. Typically, this sort of precise mapping will apply 
only when the topological weight w(C) can be expressed 
in terms of local quantities. However, once this map is 
made, the extension to more general non-unitary cases is 
obvious. 

There are two types of loop models to be discussed 
here. Fully packed loop models have the requirement that 
(at least) one loop touches every site on the lattice. Mod- 
els which do not have this requirement are referred to as 
dilute. The latter naming is somewhat cavalier: the crit- 
ical points of interest typically are believed to describe 
a phase transition between a dilute phase (where only a 
set of measure zero of the sites are touched on average by 
a loop), and a dense phase (where only a set of measure 
zero of the sites are not touched by a loop). Neverthe- 
less, since at the critical point the loop density is still 
zero (like the magnetization in the Ising model), I will 
persist with the name. 

The dilute and fully packed loop models arise here in 
very different ways. The fully packed models are found by 
expressing the transfer matrix of a lattice height model 
in terms of the generators of an algebra, and then find- 
ing a loop representation of the same algebra. It is then 
possible to define a loop model which has the identical 
partition function. The dilute models arise more indi- 
rectly: by studying the ground states of the Hamiltonian 
associated with the (same) lattice height model, one can 
infer what the domain walls are. Like in the Ising model, 
the loops are the domain walls. This argument is not rig- 
orous like that for the fully packed loops, but substantial 
consistency checks are made using exact results from the 
corner transfer matrix and the scattering matrix. 

The outline of the rest of the paper is as follows. In sec- 
tion [Til I present the main results. In section [TTTJ I review 
how the "minimal models" of conformal field theory are 
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related to both dilute and fully packed loop models. In 
section HVl I introduce height models which have critical 
points described by the superconformal minimal models 
and SU (2)2 conformal field theory. There these height 
models are used to define fully packed loop models whose 
topological weight is given by a chromatic polynomial. In 
section [V] dilute loop models are found which have these 
critical points, and then generalized to include SU(2)k 
as well. In appendix [5] the relevant conformal field theo- 
ries are reviewed briefly In appendix [B] the appropriate 
boundary conditions for the lattice models are discussed, 
while some technical details on the BMW algebra are 
collected in appendix [Cl 



II. THE LOOP MODELS AND THE RESULTS 

In this section I introduce the dilute and fully packed 
loop models studied in this paper, and describe their crit- 
ical points. The following sections are then devoted to 
explaining and deriving these results. 



A. Dilute loops 

The dilute-loop model discussed in this paper is a gen- 
eralization of the 0(n) loop model with two types of 
loops. The two types are strongly coupled: each link 
containing one kind of loop must also contain the other. 
A configuration is displayed in figure [21 where the two 
types of loops are represented by solid and dashed lines. 
The loops are drawn across the diagonals of the plaque- 
ttes of the square lattice. 




FIG. 2: A typical configuration in the new dilute loop model. 



In both this model and the 0(n) model there may be 
two loops coming into the same point on the lattice, and 
it is required that they do not cross. The new model is 



distinguished from the 0(n) model by how the loops are 
resolved when they potentially cross at the center of a 
plaquette. In the new models, it is required that loops 
of each type do not cross. Loops of different types may 
cross; the four different possibilities are illustrated below 
in figure [T5l 

Topologically-invariant quantities of such loop config- 
urations are the numbers Af and M. of each kind of loop, 
and C, the number of plaquettes with a resolved poten- 
tial crossing at their center. In figured] these are Af = 5, 
M = 4, and C — 4. The topological weight is 



w(C) 



n H m M b c 



(3) 



where n, m and b are parameters. For m = b = 1, this 
weight reduces to that of the O(n) model. It is convenient 
to parametrize n and m by 



n = 2 cos 
m = 2 cos 



7T 



p - k + 1 



k + 2 



(4) 



I will argue that with appropriate boundary conditions, 
these dilute loop models have a critical point described 
by a unitary conformal field theory when p and k are 
integers with p > k + 1 and k > 1. These conformal 
field theories are called coset models, and are described 
briefly in appendix [5] The k = 1 case reduces to the 
0(n) model, which is well known to be described by the 
conformal minimal models For k = 2, the confor- 
mal field theories are the superconformal minimal mod- 
els, whereas for p — > 00 (n — ► 2), they are the SU(2)j. 
Wess-Zumino-Witten models. It is natural to conjecture 
that this loop model has a critical point when n < 2 
and m < 2, although for p and k outside these special 
values, the corresponding conformal field theories will be 
non-unitary. 



B. Fully packed loops 

A fully packed version of the 0(n) loop model is illus- 
trated in figure [3] In this version the loops cover every 
link of the square lattice. The only degree of freedom is 
how the loops are resolved at each site of the lattice: there 
are two possible resolutions which do not allow crossings. 
The topological weight is then 



W F (C) = (y/Q) 



(5) 



when there are Af of these loops. The weight per loop 
is labeled \J~Q here because for weight per unit length 
t = 1 and appropriate boundary conditions, this model 
is equivalent to the Q-statc Potts model at its self-dual 
point. This is shown using the high-temperature ex- 
pansion of the Potts model [13, [HI], which will be re- 
viewed in section [Til Al The Potts self-dual point is crit- 
ical when Q < 4, so that the critical point occurs for a 
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FIG. 3: A typical configuration in a fully packed 0(n) loop 
model 



weight per loop < 2, just like the dilute 0(n) models. 
Also like the dilute 0(n) models, the continuum limit of 
the critical Potts models are described by the minimal 
models of conformal field theory when Q takes on spe- 
cial values. One obtains the pth minimal model when 
^/Q = 2cos(-7r/(p + 1)) with p integer; note the shift of p 
relative to the dilute value found in ((4]) with k = 1. This 
critical point is often called the dense critical point of the 
0{n) loop model 0,0. 




FIG. 4: A typical configuration in the fully packed k — 2 
model 

In terms of the more general coset conformal field the- 
ories with central charge (|A3j) . the minimal models and 
hence the dense critical points have k = 1. In section 
IIVI I derive fully packed loop models which have critical 
points corresponding to the k — 2 case, the superconfor- 
mal minimal models. A typical configuration in this fully 
packed model is illustrated in figure 2] There are three 
allowed configurations for each vertex, the two ways of 



resolving the lines so that they do not cross, and a third, 
which I call the "intersection" . Since these intersections 
are not resolved, the configurations do not really form 
loops, but for lack of a better name, I still call the de- 
grees of freedom loops. 




FIG. 5: The dual graph for the loop configuration in figure [4] 



The topological weight of the model is not defined 
by resolving the intersections to make the configurations 
into non-crossing loops. Instead, it is given in terms of 
the chromatic polynomial. To define the chromatic poly- 
nomial, it is first useful to note that the loops divide 
the plane into different regions. These regions can be 
thought of as countries on a map, and a Q-coloring gives 
each country one of Q colors, such that adjacent coun- 
tries must have different colors. If two regions touch only 
at a point (or any countable set of points) , they are not 
adjacent. This is convenient to formulate in terms of the 
dual graph. Each region corresponds to a vertex on the 
dual graph, and two vertices are joined by an edge if the 
corresponding regions are adjacent. The dual graph for 
the loop configuration in figure 0] is illustrated in [51 Note 
that the dual graph only depends on the topology of the 
original loop configuration. The number of Q-colorings 
Xq (£) of the dual graph is the number of ways of assign- 
ing Q colors to each vertex with the constraint that any 
vertices joined by an edge must be a different color. 

There exists no closed-form formula for the number 
of colorings of a general graph. However, \Q can be 
determined recursively. Consider two vertices connected 
by an edge I (i.e. two regions sharing a boundary in the 
original picture). Then define T>iC to be the graph with 
the I edge deleted, and CiC to be the graph with the two 
vertices connected by I joined into one. Then it is simple 
to prove that 

Xq {£)=Xq{'D 1 L)-Xq{C 1 C). (6) 

In the dual graph of C the two vertices connected by I 
must be colored differently, while in T>iC, they may be 
the same, Thus x<g(X>z£) > xq(£)> an d the number of 
colorings overcounted is XQ (£/£)■ The number of col- 
orings can then be determined by using this recursion 
relation to remove edges and vertices from C until it is 
a sum of terms which have no edges at all. A term with 
V vertices and no edges has \Q ~ Q V ■ An illuminating 
example is the case where there are no intersections, so 
that the configurations really are loops. The dual graph 
is then tree-like (it has no faces). Vertices at the ends of 
the tree can be removed one by one using ©, with each 
vertex being removed giving a factor Q — 1. This yields 
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Xq = Q(Q — l) v_1 for a tree with V vertices. Using the 
recursion relation ([B]) gives for the loop configuration in 
figure 0] and the dual graph in figure [5] 

x Q = Q(Q-i) 9 (Q 2 -3Q + 3) 2 . 

Repeatedly applying the recursion relation gives xq (£) 
as a polynomial in Q of order V. This is known as 
the chromatic polynomial, and depends only on the dual 
graph of C. This polynomial provides a well-defined ex- 
tension of xq away from Q integer. The fully packed 
loop model is therefore uniquely defined for any Q with 
topological weight 

Afx is the number of intersections; the factor involving it 
arises from the detailed analysis below. The full partition 
function is given by using this topological weight with 
([T8]) . where the weight per unit length of loop is t = 1. 
The sum is over all loops on the square lattice such that 
each plaquette has one of the three configurations drawn 
in figure [T2l 

In section HVBl I prove that the loop model with topo- 
logical weight ([7|) has a critical point described by the 
pth superconformal minimal model (i.e. those in (|A2|) 
with k = 2), when 

Q = w(-^) . (8) 

These fully packed models presumably remain critical for 
all real p, but are described by a non-unitary conformal 
field theory for p non-integer. 

This loop model is very similar to that of [121 ]. which 
arises from the low-temperature expansion of the Q-state 
Potts model. The topological weight is given by the chro- 
matic polynomial. However, the loops there are dilute, 
not fully packed, and intersections with three lines com- 
ing out are allowed there. As a result, the critical point 
there is described by the conformal minimal models in- 
stead of the superconformal minimal models. 

This is not the only fully packed loop model which 
has critical points corresponding to the superconformal 
minimal models. One based on non-intersecting loops on 
the "copper-oxide" lattice will also be briefly discussed 
below. 



III. THE O(n) LOOP MODELS 

The strategy of this paper is to use known integrable 
lattice models to find loop models with critical points 
described by conformal field theories. In this section this 
strategy is used to connect the 0(n) loop models with 
the minimal models of conformal field theory. None of 
the results in this section are new, but possibly some of 
the arguments are. After introducing the lattice models, 
the two different possible loop expansions are derived. 



Even though both loop expansions come from the same 
lattice model, they arise in very different ways. The fully 
packed loops arise by rewriting the model in terms of 
an algebraic formulation of the transfer matrix, which in 
most contexts is a high-temperature expansion. Dilute 
loops arise from a domain- wall expansion, which is a low- 
temperature expansion. 



A. RSOS lattice models 

The lattice models here are the "restricted-solid-on- 
solid" (RSOS) models introduced by Andrews, Baxter 
and Forrester [l3j]. The degrees of freedom are integer- 
valued heights hi on each site i of the square lattice. 
There are two types of restriction. The first is that 
hi only takes a finite number of possibilities, which are 
numbered from 1 . . . p for some integer p. The second is 
that heights on adjacent sites must differ only by 1, i.e. 
\hi — hj\ = 1 for i and j nearest neighbors. The latter 
restriction will be modified for the models in later sec- 
tions. In order to utilize crucial exact results, I study only 
height and loop models on the square lattice. However, 
frequently (but not always), one obtains similar results 
on different two-dimensional lattices. 

The p = 3 case is the Ising model: on one sublattice the 
height must always be 2, so that on the other sublattice 
the heights 1 and 3 play the roles of the + and — Ising 
spins. The p = 4 RSOS model is known as the hard 
square model [11] : a typical configuration is illustrated 
in figure [5] below. 

The RSOS models are interesting and useful because 
for special choices of the Boltzmann weights, they are 
integrable. Integrability allows a number of important 
physical quantities to be computed exactly, under a set 
of standard analyticity assumptions. Critical points and 
the associated exponents can be found using a powerful 
technique called the corner transfer matrix [ll[ . Confor- 
mal field theory also provides a list of critical exponents 
associated with critical points. By comparing the two 
lists, one can usually identify the conformal field theory 
describing the continuum limit of any critical point in 
an RSOS model. For these RSOS models, there are two 
known critical points for each value of p > 3. This pa- 
per is concerned with just one of them, the one separating 
regimes III and IV in the nomenclature of [l3| . This crit- 
ical point is described by the »th minimal model, with 
central charge given by (|A1[) [l6j| . 

Finding the loop models associated with this RSOS 
model and hence the conformal minimal models requires 
analyzing the Boltzmann weights in detail. The weights 
contains nearest-neighbor and next-nearest-neighbor in- 
teractions, so the total weight can be written as a prod- 
uct of weights assigned to the four heights around each 
square plaquette. Such a model is called "interaction 
round a face" [ll|. The transfer matrix T^ji' is defined 
to act across the diagonals of the lattice, as illustrated in 
figure [HI It therefore acts on a zig-zag row of heights hi, 
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FIG. 6: How the transfer matrix acts 



properties of the model follow solely from the algebra, 
and are independent of the representation of the ej. 
With appropriately-chosen boundary conditions, the par- 
tition function is independent of representation. This is 
shown for a two-dimensional surface with the topology 
of a sphere in the appendix [Bj while detailed discussions 
for the cylinder and torus can be found in [l?], [Hj]. In 
many considerations, understanding the boundary condi- 
tions precisely is very important; for example, the central 
charge of the continuum conformal field theory depends 
on them. For the results here, the details of the bound- 
ary conditions are not particularly important - what is 
of main importance is that appropriate ones exist. 



i = .. . 2N, taking it to the next row, labeled by h^. It 
can be broken into pieces %, which depend only on the 
heights around a square, and can change hi to h[. Then 
Th h' can be written in the form 



T = %T 3 



Tin -1I0T.T 



where the identity matrix Ii = 5hih' 
critical point of interest here, we have 



T = L 



T2N-2I2N (9) 
. For the RSOS 

(10) 



where x is a parameter, and the elements of the matrix 
ei for i even are [TH, [TH 



= 5, 



hi-ih i+ i " 



ii'i 



i+l\q 



n 



(ii) 



where 



h} q = ( q »- q -»)/(q- q -i), 



3 even, &i 



Acting with ej therefore allows the height on the ith site 
to change when hi-i — hi+i. For i odd, the matrix e, is 
essentially the same: 



»J?lAi]g 



n 



(12) 



odd, jV« 



The RSOS transfer matrix written in this form exhibits 
an interesting and important algebraic structure: the 
generate the Temperley-Lieb algebra [Hj]. It is straight- 
forward to check that they satisfy 



e? = (q + i l ) e u 

e l ej = ej e 4 (\j — i\> 2). 



(13) 



This algebra first arose in studies of the Q-state Potts 
model, whose transfer matrix also can be written in the 
form (|9ll0p . with the satisfying the same algebra p^]) 
with q + q^ 1 = \[Q. Since then, it has been shown how 
to write many other lattice models in this same form. 

Writing a transfer matrix in terms of the Temperley- 
Lieb generators is exceptionally useful because many 



B. Fully packed loops 

This construction of the fully packed loop model asso- 
ciated with each critical RSOS model is virtually identi- 
cal to that of the fully packed loop model for the Potts 
model [13, EH • This is easiest to do by using a graphical 
representation of the Temperley-Lieb generators [17[ • 

Fully packed loop configurations are not in one-to-one 
correspondence with those of the height or Potts model. 
Rather, for each height configuration, the product in (J9j) 
is expanded out by using (jTU|) . One obtains either the 
identity / or the Temperley-Lieb generator e for each 
plaquette, so that for V plaquettes there are 2 V terms. 
Each term can be graphically represented by drawing one 
of the two pictures in figure [7] for each plaquette, de- 
pending on whether I or e appeared for this plaquette. 
This gives a fully packed configuration of self- avoiding 
and mutually-avoiding loops on the dual square lattice. 
Every link of the dual lattice is covered by a loop, and 
at each vertex there are two possible ways of resolving 
the lines so that the loops do not cross. These properties 
can be preserved at the boundaries as well, as described 
in appendix IBI 




FIG. 7: The graphical representation of the Temperley-Lieb 
generators 



To obtain the fully packed loop model, one must sum 
over the heights in the RSOS representation (or the spins 
in the Potts representation) , leaving only the loops them- 
selves as the degrees of freedom. The sum is not diffi- 
cult to do directly in the RSOS representation (or 
in the Potts representation [l(| Hl|). by using the ex- 
plicit representation for e, as given in pip . It is more 
instructive, however, to not utilize the details of any ex- 
plicit representation, but rather to exploit the properties 
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of the Temperley-Lieb algebra. In the transfer matrix, 
summing over heights amounts to multiplying the ma- 
trices defined by Ii and e,. Therefore, the algebra itself 
can be represented graphically as in figure [5] From this 
it is apparent that small closed loops receive a weight 
q + q^ 1 . By using both relations, it is easy to check 
that all loops receive this same weight q + q^ 1 , no mat- 
ter what their size (for details, see appendix |B|) . Setting 
x = 1 to make the model isotropic (invariant under 90- 
degree rotations) yields the fully packed 0(n) loop model 
discussed in the introduction. The weight per loop is 
n = q + q^ 1 = 2cos[7r/(p + 1)] for the RSOS model, or 
n = y/Q for the Q-state Potts model. The weight per 
unit length of loop t in this case is unimportant, because 
every link is covered by a loop. 



xx 




xx 




xx 
xx 



FIG. 8: The graphical representation of the Temperley-Lieb 
algebra 



The 0(n) loop model can of course be defined for any 
n, but is only critical for \n\ < 2. This is shown by using 
the representation of yet another lattice model, the six- 
vertex model, in terms of the Temperley-Lieb algebra 
[Til ] . Another interesting thing to note is that even when 
| ?i | < 2, the conformal field theory is unitary only for 
integer p > 3. The representation pip exists for any 
value of q — e i7r /(p +1 ), but the only ones leading to real 
and positive Boltzmann weights for q + q^ 1 < 2 occur at 
p integer. 



C. Domain walls in the RSOS model 

The dilute 0(n) loop model also has a critical point 
related to the conformal minimal models Q. A way of 
showing this heuristically is to compare the Coulomb-gas 
description of the minimal models [5|, to that arising 
from the lattice model [l| . A way of showing this directly 
is to study a different series of integrable RSOS models 
known as the "dilute Temperley-Lieb" models [2(| HH, 
[22| . In this section I give another method of relating the 
two, by using the RSOS models defined in section Ull Al 
This is the method which readily generalizes to k > 1. 

It is useful to start by discussing how the dilute 
Temperley-Lieb models are related to the 0(n) loop 
model [20j, |2l| . The dilute models are defined in terms of 
integer heights on the square lattice like the RSOS mod- 
els discussed above in section IIII Al but the restriction 
on neighboring heights in the dilute models is relaxed 



refers to the fact that because hi — hj can vanish for near- 
est neighbors, there need not be domain walls on every 
link. Just drawing domain walls does not immediately 
make it into the 0(n) loop model: there is no reason a 
priori why a region of heights 1 and 2 has the same weight 
as a region of heights 2 and 3. Nevertheless, when the 
Boltzmann weights are tuned appropriately, one can per- 
form the sum over heights to leave a model consisting 
of self-avoiding and mutually-avoiding loops with weight 
n per loop [20j. By then doing a corner transfer matrix 
computation, one finds indeed that the continuum limit is 
described by a conformal minimal model, with the weight 
n per loop related to p by n = 2cos(tt/p) (note the shift 
in p as compared to the fully packed case) [13] • 

The loops discussed in the remainder of this section 
can be thought of as domain walls for the RSOS model 
defined in section lTlI Al However, they are not simply do- 
main walls between differing heights, as they are in the 
dilute model. Instead, they are best thought of as domain 
walls between "ground states" in a 1+1-dimensional pic- 
ture. Namely, one chooses a Euclidean time direction, 
and takes the continuum limit in the time direction to 
obtain a one-dimensional quantum Hamiltonian, and find 
its ground states. Excitations in the one-dimensional pic- 
ture are then kinks, or defects, which separate regions 
comprised of the different ground states. The world lines 
of the kinks in 1+1 dimensions can be thought of as do- 
main walls in the two-dimensional classical model. Ob- 
viously, this picture is heuristic, but below I will explain 
how to use it to develop exact results for the RSOS model 
into a precise conjecture for a dilute loop model. 

A quantum Hamiltonian for the RSOS models is ob- 
tained in the limit x — * 0. One finds that the critical 
properties of the model are independent of along as 
its sign is not changed. Thus taking x — > in ^ and 
(I10[) gives a quantum Hamiltonian acting in the same di- 
rection as the transfer matrix, i.e. across the diagonals of 
the lattice. It is (for x > 0) 



2N-1 



(14) 



to hi 



0,±1. The loops are simply the domain 



walls between different heights. The "dilute" in the name 



This Hamiltonian acts on the "height chain" 
(ho, h-2,..., h2N)- 

Using (fTTj) , it is easy to find candidates for the ground 
states of (fl4l . The matrix elements of are positive 
when hi-i = hj+i, and vanish when =/= ftj+i- Thus it 
is natural to expect that a ground state will be dominated 
by states that obey hi+\ = hi-i for all i. The RSOS 
models defined in section IIII Al have the restriction that 
hi — hi-i = ±1 for all i. Thus potential ground states 
here are of the form (. . . r, r + 1, r, r + 1, . . . ). 

Finding which one or ones of these potential ground 
states dominates the actual ground state or states is in 
general an imposing problem. For the integrable Hamil- 
tonian (|14p. however, it is possible to answer this question 
by using the corner transfer matrix technique [111 [l3l ]. 
The result is that there is a ground state dominated by 
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each of the potential ground states. Thus this critical 
point is a multicritical point, where all the potential or- 
derings coexist fl6j |. 

I illustrate this here in the simplest cases, and then in 
the next subsection IIII Dl explain how the general result 
follows from the corner transfer matrix. The simplest 
case is the Ising model, which has p = 3. The restriction 
that \h% — hi+i\ = 1 means that every other site on the 
chain (say the even sites) must be occupied by height 2. 
The odd sites are occupied by heights 1 and 3, which 
play the role of the spins ± in the usual formulation of 
the Ising chain. Using the explicit matrix elements given 
in (fTTj) means that — a for even i = 2j yields a potential 
energy which is — v2 if /laj-i = ^2j+i, and for /i2j_i 
/i2j+i- For odd i = 2j — 1, 



e 2 j-i 



J_/l 1 
V2 V 1 



where the rows and columns are indexed by h' 2 j_ 1 = 1,3 
and/i2j-i = Irrespectively. Rewriting the Hamiltonian 
in terms of the Pauli matrices <J a (j) gives 



N/2 



where t = 1. This is exactly the ferromagnetic Ising 
spin chain; t — 1 is the critical point. For t < 1, the 
chain is ordered: there are two ground states, one domi- 
nated by the configuration with by all spins up, the other 
dominated by the configuration with all spins down. In 
the height language, these correspond respectively to the 
configurations (121212 . . . ) and (323232 . . . ). For t > 1, 
the model is disordered. Thus t = 1 is the well-known 
Ising critical point between the ordered and disordered 
phases. Excited states in the ordered phase contain de- 
fects, i.e. states which contain regions of both ground 
states locally. The equivalence of the Ising model to the 
dilute loop model with n = 2(cos(-7r/3)) = 1 was already 
discussed in the introduction. In the height language, 
the loops are domain walls separating regions of height 1 
and height 3. 

The two ground states in the p — 3 case are symmetry- 
equivalent: the critical point is just the usual one where 
the discrete symmetry becomes no longer spontaneously 
broken. For p > 3, the different ground states are no 
longer all related by any obvious symmetry, but the fine 
tuning necessary to get to the integrable multi-critical 
point with Hamiltonian (|14p makes them all degenerate. 

To see this in more detail, consider the p = 4 case, the 
"hard-square" model [ill ]. The hard-square tiles have 
diagonal length twice the lattice spacing, and are placed 
with their centers on the sites of a square lattice, with 
the rule that no two tiles can overlap. This then forbids 
tiles from being adjacent to each other. In the height 
language, this means that heights 1 and 4 correspond 
to the tiles, and 2 and 3 correspond to empty sites (on 
one sublattice the height is always even, on the other it 




FIG. 9: A typical configuration in the RSOS height model 
with p = 4. The lines are the domain walls corresponding to 
dilute loops, drawn so that they do not cover the heights that 
they go through. 



is always odd, and the Boltzmann weights are invariant 
under h — > p+l — h.) The restriction that \hi — hj\ = 1 for 
i next to j enforces the hard-square restriction. A typical 
configuration is displayed in figured The Hamiltonian is 
easiest to write by treating the square tile as a hard-core 
boson created by d\ . It is then [23[ 

H = [~W (ck + dt) (1 - ni_i)(l - n l+1 ) 

i 

+U ni + V ni n i+2 ] (15) 

where rii = d\di is the number operator at site i. The 
critical (fT4|) corresponds to the values 



-1/2 



U c = 2 + l/<f> 2 (16) 



where = 2 cos(7r/5) = (1 + y/E)/2 is the golden mean. 

To understand the multiple ground states, it is useful 
to first take the limit where |f/W| and IU/I4 7 ! are large, 
so that the off-diagonal term proportional to W can be 
neglected. For U > — V then the ground state is the 
state with no bosons, while for U < — V there are two 
ground states, each of which has a particle on every other 
site. In the latter regime, a Z2 symmetry (translation by 
one site) is spontaneously broken. At U = — V, there 
is a first-order phase transition between the Z2-broken 
phase to an unbroken phase. In height language the 
two ground states in the Z 2 -broken phase are dominated 
by the states (. . . 12121212 . . . ) and (. . . 34343434 ...) re- 
spectively, while the single ground state in the unbroken 
phase (the empty state in the bosonic language) is dom- 
inated by (. . . 32323232 . . . ). This first-order transition 
persists as W is included; doing perturbation theory gives 
the location of this transition at next order to be [231 



U 
V 



V 2 '' 



(17) 



9 



Although the perturbative computation is only reliable 
to order (W/V) 2 , it turns out that (fTT|) is the exact first- 
order transition line for U > U c . Along this line the 
two phases coexist - the hard-square model has three 
symmetry-equivalent ground states. When the couplings 
are tuned to the values in (fT6|h this first-order line ter- 
minates in the tricritical point [111, IH| ■ 

This information about the ground states is sufficient 
to guess what the lowest-energy excitations are; in the 
next subsection IIII Dl I will justify these guesses us- 
ing the corner transfer matrix. Excited states along 
this first-order transition line are defects, where two 
different ground states meet. For example, the state 
(. . . 32323232121212323232 . . .) has two defects, and en- 
ergy order U higher than the ground state. Thus each 
of these defects has energy of order U/2. Such a de- 
fect occurs at site i any time — = 2. There 
are two kinds of fundamental defects, with the same en- 
ergy. Denoting the three ground states as G12, G23 and 
G34, where G a b — (• • • abababab . . . ) , one kind of fun- 
damental defect separates G12 and G23, while the other 
separates G23 and G34. The former defect is located at 
the site where there is a height 2 between the heights 
1 and 3, while the latter is located at the site where 
there is a height 3 between heights 2 and 4. A defect 
(. . . 1212343434 . . . ) separating G12 from G34 has twice 
the energy as the fundamental defects, and is comprised 
of two fundamental defects. 

These statements can be translated from the one- 
dimensional quantum Hamiltonian back to the two- 
dimensional hard-square model (24|. The three ground 
states of the Hamiltonian result in three different ways 
the free energy density can be minimized. One can easily 
check that at the critical point, the weight of the config- 
uration with all heights 1 and 2 (which in a slight abuse 
of notation can be labeled G12) is larger than that of 
the configuration containing all heights 2 and 3 (the 2d 
analog of G23). On the first-order line, there are three 
different ways the free-energy density can be minimized. 
The reason is that there are more configurations "close" 
to G23 than there are to G12 or G34, so the increased 
entropy will compensate for the increased energy. This 
is the cause as well of the shift of the first-order line in 
(fT7|) . With G23, one can change any height 2 to height 
4, and any height 3 to height 1, without violating the 
hard-square rule. With G12, one can change any height 

1 to height 3, but it is not possible to change any heights 

2 to height 4 without violating the restriction, unless one 
also changes the four heights 1 around it. In the lan- 
guage of hard squares, there are more ways of adding a 
square to the empty configuration than there are of re- 
moving one from the full configuration, since in the latter 
the squares occupy only every other site. Thus the three 
ground states turn into the three minima of the free en- 
ergy. 

The partition function can therefore be rewritten as 
a sum over domain-wall configurations separating the 
free-energy minima. The excited states of the Hamilto- 



nian become the domain walls in the 2d classical model, 
separating regions of Z2 broken and unbroken symme- 
try. These domain walls are the world-lines of the one- 
dimensional defects, and their contribution to the free 
energy depends only on their length times some positive 
constant fr>. By the symmetry h — > p+1 — h, fu must be 
the same for both kinds of domain wall. A non- vanishing 
fo is the reason the transition is first-order: the tricrit- 
ical point occurs when the couplings are tuned to make 
fo — > 0. In this limit there is no free-energy penalty for 
creating a domain wall, and they proliferate. 

It is finally possible to give the argument as to why the 
hard-square model is equivalent to a loop model, in the 
sense that they are described by the same field theories in 
the continuum limit. I have argued that the hard-square 
model along its first-order line can be thought effectively 
as having three free-energy minima. Domain walls sepa- 
rating regions of these minima have a free-energy weight 
per unit length. The domain walls form loops which do 
not touch: as noted above, domain walls separating G12 
and G23 are located along sites with height 2, while do- 
main walls separating regions G23 and G34 are located 
along sites with height 3. These domain walls are illus- 
trated in figure I emphasize that this picture involving 
free energies and domain walls is effective: these should 
be thought of as renormalized quantities. 

Nevertheless, these results imply that along the first- 
order line, the partition function of the hard-square 
model can be written in the form 

Z p=4 « 5>,(/;)e- /DL . (18) 
c 

£ is a configuration of loops of total length L on the links 
of the dual square lattice, with the condition that they 
cannot touch. In this expansion, no distinction is made 
between the two different types of domain walls: the fact 
there are two types of loops can be taken into account in 
the weight w(C). To do this, draw a loop configuration. 
The restrictions on heights mean that a loop separates a 
region of G23 from a region of either G12 or G34 . Since the 
energies of G34 and G12 are related by the Z2 symmetry, 
we can then sum over these two possibilities. Half the 
regions inside the loops therefore must be G23, while the 
other half can be either G12 or G34. This results in a 
weight w(C) = 2 jV ^ 2 , where Af is the total number of 
loops, so that the weight per loop is y/2. This is the 
dilute 0(n) loop model with n = The (tri)critical 
point occurs at fjj = where there is no penalty for 
longer loops and they can proliferate. 

The loop model therefore arises here indirectly, as op- 
posed to the dilute Temperley-Lieb models, where the 
domain walls and hence the loops are manifest in the 
original lattice height configurations [2(| Hl[ . All argu- 
ments give the same result: the conformal field theory 
with p = 4 describes the critical point of the 0(^/2) di- 
lute loop model on the sphere. 
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D. Dilute loops and the corner transfer matrix 



In the preceding subsection IIII CI I detailed an argu- 
ment relating the RSOS model with p = 4 to the 0{-\/2) 
loop model. The argument can be strengthened and gen- 
eralized to arbitrary p by using the corner transfer matrix 
technique [ill [l3| . The corner transfer matrix allows one 
to find the ground states in general, and also gives essen- 
tial information about the defect energies. 



'b 



m h 



b c 



FIG. 10: The corner transfer matrix computes the partition 
function in one quadrant of the square. In the frozen limit, the 
spins are all the same along a given diagonal of each quadrant. 



Consider a two-dimensional lattice model on a square 
lattice, as illustrated in figure [TO] The corner transfer 
matrix acts on TV sites along a line, taking them to N 
sites along a line perpendicular to the original line. The 
partition function for the full model is then built up from 
four corner transfer matrices. 

The glory of the corner transfer matrix is that, as op- 
posed to the usual transfer matrix, its eigenvalues can 
be determined exactly and explicitly if the model is inte- 
grable 11| . Its eigenvalues have special analyticity prop- 
erties in one of the couplings (x here) , which allows one 
to show that they must be independent of other cou- 
plings. In the hard-square model, the eigenvalues of the 
corner transfer matrix are the same along the entire first- 
order line. Therefore, they can be computed by going 
to the far "end" of this line, where in the Hamiltonian 
\U/w\ — \V/w\ — > oo. In this limit, not only are the 
three states G12, G23 and G34 the exact ground states 
of the Hamiltonian, but every height configuration is an 
eigenstate. This limit is "frozen" , because in the two- 
dimensional lattice model all the heights along a diagonal 
are frozen to be the same. For example, in the northeast 



and southwest quadrants of figure [101 the heights along 
each diagonal in the northwest/southeast direction must 
be the same. Moving away from the frozen limit, the 
eigenstates of the corner transfer matrix will change, but 
the eigenvalues, and hence the type of ground states, stay 
the same all along the first-order line. The same picture 
holds for all p: there is a first-order transition line going 
from a multicritical point to a frozen point. 

From the corner transfer matrix eigenvalues one can 
compute the local height probabilities P(a\b,c), defined 
as the probability that the height in the center will be 
height a, given that the heights around the edges of the 
square are b and c, as illustrated in figure [TOl The N 
heights along a line starting at the center are labeled hi, 
so that hi = a, /i/v-i = b, and h^ — c. Then 13| 



P(a\b,c) =V a 



{h 2 ,h 3 ,...hf 



N-2 



E(h 1: h 2 ,...,h N ) 



E{h 1 ,h 2 , . . -,h N ) = ^ 3\ h 3 ~ h 3+z\ 



(19) 



(20) 



where V a is a (known) factor independent of the heights 
other than hi = a, and q parametrizes the first-order 
line, so that q — > is the frozen limit and q — > 1 is the 
multicritical point. The sum is over all allowed height 
configurations, and the contribution of each distinct con- 
figuration to this probability is obvious. The j in the 
expression for E arises because there are j sites along 
the jth diagonal, and in the frozen limit each site con- 
tributes equally to the eigenvalue. 

The results found in the previous subsection llll Cl using 
the Hamiltonian are therefore identical to those implied 
by the corner transfer matrix, and generalize to all p. 
Ground states in the frozen limit have E = in (|20[) . 
These are states where the heights on every other site 
are the same, and it is natural to expect that the ground 
states all along the first-order line will be dominated by 
these configurations. Because of the restriction that \hj — 
hj + i\ = 1, all states of the form Gjj+i = + 
1) jj 3 + 1 • • • )■ F° r the fth model, this means there are 
p — 1 ground states. A defect at site j + 1 occurs when 
\hj — hj + 2\ = 2; such a defect runs all along the diagonal, 
and so has an energy of 2 per unit length in the units 
implied by (|19I20|) . All defects have this energy per unit 
length, no matter which ground states they separate. 

The translation of these results to the two-dimensional 
lattice model is essentially the same as for p = 4 [lj|. 
Each of the p — 1 ground states Gjj+i corresponds to a 
degenerate minimum of the free energy. Thus for q — ► 
the partition function is dominated by regions of these 
minima, separated by domain walls which all have the 
same energy per unit length. For general p, as opposed 
to p = 3,4, there is no symmetry forcing the domain walls 
to have the same energy, but rather it is a consequence of 
the fine tunings necessary to make the model integrable 
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and on the first-order transition line. As q — > 1, the mul- 
ticritical point of interest is approached, and the energy 
per unit length of the domain walls goes to zero. 

The restrictions on adjacent heights and on domain 
walls can be simply encoded in terms of adjacency di- 
agrams. For heights, the diagram has a node for every 
allowed height, and two nodes are connected by a edge 
if the corresponding heights are allowed to be on adja- 
cent sites. For ground states, the diagram has a node 
for every ground state, and two nodes are connected if 
a fundamental domain wall can separate the two ground 
states. The RSOS models discussed above have p — 1 
different ground states G12, G23, G34,.... The funda- 
mental domain walls separate successive ground states 
on this list: i.e. a region of Gh,h+i can only be adjacent 
to Gh-i.h and Gh+i,h+2- Thus in this case the adjacency 
diagrams for the heights and for the ground states look 
the same, except there are p nodes in the height dia- 
gram, and p — 1 nodes in the ground-state diagram. For 
both types, the jth node is connected to the j + 1 and 
j — 1 nodes. The nodes of the ground-state diagram cor- 
respond to the edges of the height diagram; this remains 
true for the more general RSOS models discussed in the 
next section. The two diagrams for p = 5 are displayed 
in figure [TT1 



G 12 G 23 G 34 G 45 



FIG. 11: On the left, the adjacency diagram for the heights of 
the p — 5 RSOS model; on the right, the adjacency diagram 
for the ground states of the same model. 



By the same argument as for p — 4, these results make 
it likely that the pth minimal conformal field theory de- 
scribes the continuum limit of the critical point of a loop 
model with partition function of the form pil8[) : 



w(C)e 



-fnL 



(21) 



with suitable boundary conditions. To obtain the topo- 
logical weight w(C), one must sum over the all the height 
configurations consistent with a fixed loop configuration. 

The simplest way to do this sum is to define the ad- 
jacency matrix A for the ground states. This is the ad- 
jacency diagram in matrix form: each row and column 
of this matrix corresponds to a ground state, with the 
A a b = 1 if the ground states labeled by a and b can be 
separated by a fundamental domain wall, and zero oth- 
erwise. For the RSOS models here, it is 

Ars &\r— s|,l (22) 
for r, s = 1, . . . ,p — 1. For example, for p = 5 it is 



A = 



/0 1 0\ 
10 10 
10 1 

v0 1 0/ 



The topological weight w(C) is the number of height con- 
figurations possible for each loop configuration C. This 
is easily written in terms of A. Index each loop by j, 
such that it separates regions of ground states labeled 
by a(j) and b(j). Then the number of allowed height 
configurations is 

where the sum is over all ground states in each region. 
Thus adding another loop amounts to multiplying by 
A. For a large number of non-intersecting loops TV, one 
therefore obtains 

w(C) ra 

where A is the largest eigenvalue of A. For the adjacency 
matrix (l22l). this is 



A = 2 cos 



Thus the RSOS models in this section indeed correspond 
in the continuum limit to the 0(n) loop model with n = 
2cos(ir/p). This agrees with the two earlier cases we 
derived: n = 1 for the Ising model, and n = \/2 for the 
hard-square model. 

A very nice corroboration of this picture comes from 
studying the scattering matrix of the 1+1 dimensional 
theory in the continuum limit of the first-order line [25| . 
In field-theory language, moving along the first-order line 
away from the multi-critical point corresponds to a per- 
turbation by the relevant $13 "energy" operator. This 
field theory is integrable like the underlying lattice model, 
so one find the quasiparticles and their exact scattering 
matrix. These quasiparticles turn out to be kinks which 
are defects between p — 1 ground states, exactly as we 
have seen. This identification of the worldlines of these 
kinks as the loops in the 0(n) model with n = 2 cos(7r/p) 
was made long ago 26]. In fact, the correspondence is 
even closer: the scattering matrix itself can be expressed 
in terms of the Temperley-Lieb algebra 27]. Namely, 
there is a representation of / and e acting on the kinks, 
such that the two-particle scattering matrix for two par- 
ticles of momentum p\ and pi is 



S{p\,P2) = f(pi,P2) (I + ff(pi,p 2 )e) 



(23) 



where / and g are known functions. This correspondence 
between the loops and the world-lines of the kinks in the 
1+1-dimcnsional description was exploited in [121 ] to dis- 
cuss different loop models for the Q-state Potts models. 



IV. NEW FULLY PACKED LOOP MODELS 

In section Mil I gave a dilute and fully packed loop 
model associated with each RSOS height model of An- 
drews, Baxter and Forrester. The conformal minimal 
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models describe the critical points. In this and the next 
section the methods developed in the last section are used 
to find new loop models from more complicated RSOS 
height models. It is then natural to conjecture that the 
conformal field theories describing the critical points of 
these height models also describe critical points in these 
loop models. 

A. Fused RSOS models 

The Boltzmann weights of the integrable lattice models 
studied here satisfy the Yang-Baxter equation. This very 
strong constraint allows many properties of the models 
to be computed exactly, like the spectrum of the corner 
transfer matrix discussed above. 

Another thing the Yang-Baxter equation allows is a 
way to construct new integrable lattice models from 
known ones. This procedure is called fusion, and was 
invented in (28|, and applied to height models in [29| . 
The fused models obtained from the RSOS models are 
labeled by an integer k, and each model within a series 
is labeled by the same integer p as before. Thanks to 
corner-transfer matrix computations 30], the conformal 
field theories describing the critical points of the fused 
RSOS models are known: they are the coset conformal 
field theories (|A2[) . which have central charge (|A3[) . 

The states of the general integrable RSOS models 
are easiest to understand in terms of representations of 
a quantum-group algebra. A quantum-group algebra 
U q (G) is a one-parameter deformation of a simple Lie al- 
gebra G. For generic values of the parameter q, U q (sl(2)) 
has irreducible spin-s representations corresponding to 
those of ordinary sZ(2), but when q is a root of unity 
such that g 2 (p +1 ) = 1, only those with s < p/2 are irre- 
ducible. The tensor product of the spin-s representation 
with the spin-1/2 one is 

(s) <g> (1/2) = (s + 1/2) + (s - 1/2) (24) 

as long as all representations involved have s < p/2, so 
e.g. {(p - l)/2) (1/2) = (p/2 - 1). Except for this 
truncation, this is the same rule as for ordinary sl(2). 

Each allowed height h of a U q (sl(2)) RSOS model 
corresponds to the irreducible representation of spin 
(h — l)/2, so the heights run from 1, . . .p. Each fused 
model is labeled by an integer k. In the fcth model, 
the height hi is allowed to be adjacent to a height h 2 if 
the corresponding representations of U q (sl(2)) have the 
tensor-product decomposition 




This can be interpreted physically by saying that the link 
variables of the fused models have spin k/2, with the orig- 
inal RSOS models discussed in section lTlI Al having k = 1. 
For k — 1 , using the tensor product (|2"4")l indeed gives the 
adjacency diagram illustrated in figure 111! The tensor 



products of higher-spin representations of U q (sl(2)) can 
be built up from ([H]). Tne fused RSOS models to be dis- 
cussed in this section have k — 2, so the links have spin 
1. The tensor product of a spin-1 representation can be 
found by using the fact that (1/2) <g> (1/2) = (0) + (1). 
Thus 

(s) <8 (1) = 00 <8> (1/2) <8 (1/2) - (s). 

For s < p/2 — 1, the tensor product of (s) with (1) is thus 
the same as for ordinary sl(2). The different rules come 
from the restriction that s < p/2. One has (p/2 — 1) <E> 
(1) = (p/2 - 1) + (p/2 - 2), while ((p - l)/2) ® (1) = 
((p — 3)/2). Note that ((p — l)/2) does not appear on 
the right side of the latter, despite its being an allowed 
representation. 

For the general model labeled by k and p, continuing 
in this fashion gives the constraints on nearest-neighbor 
heights hi and hj to be [29| 

hi = 1,2. ..p 

hi - hj = —k, —k + 2,...,k (26) 
k + 1 < hi + hj < 2p - k + 1 

These rules are symmetric under the interchange h — > 
p + 1 — h. The Boltzmann weights for k = 1 satisfy this 
symmetry as well, so those for higher k found by fusion 
satisfy this as well. For k even, all heights in a given 
configuration must be even or odd. It then is consistent 
to restrict them all to be even, because one obtains the 
same field theory for even or odd. 

B. Fully packed loop models for k = 2 

I show here how to derive fully packed loop models for 
the k — 2 fused RSOS models. This is done in a similar 
fashion as for the k = 1 case, by rewriting the transfer 
matrix in terms of the generators of an algebra, and then 
finding a graphical representation generalizing that of the 
Temperley-Lieb algebra. 

The transfer matrix for all k can be written in the form 
©. For k = 2, the weights for each plaquette are then 

% = h+ xXi + yE % (27) 

where x and y are parameters related by 



y/Q = q + q- 1 = 2 co S (n/(p + 1)) . 

The model is isotropic when x = 1/(VQ + 1) and y = 1. 
The Hamiltonian 

H = -^2(Xi + Ei) (28) 

i 

follows from the x — > limit. The representation of Xi 
and Ei in terms of the heights is given in appendix [C] 



There is one important difference between Xi and the 
Temperley-Lieb generators: can be non-zero when 
both pairs of heights across a plaquette are different (i.e. 
hi 7^ h[ and hi-\ ^ hi+x). This has important conse- 
quences for both the dilute and fully packed loop models. 

The lattice model with weights (f27|) is integrable when 
the Xi and Ei satisfy an algebra known as the 5*0(3) 
Birman-Murakami-Wenzl (BMW) algebra This al- 
gebra is best thought of as the spin-1 generalization of 
the Temperley-Lieb algebra. In fact, any representation 
of the Temperley-Lieb algebra yields one of the 5*0(3) 
BMW algebra [32] ; expressions for Xi and Ei in terms of 
the ei are given in appendix [UJ The algebraic relations 
of the Xi and E t can be found in [32j], and are given in 
appendix [Cl 

As with lattice models based on the Temperley-Lieb 
algebra, many of the properties of the models with Boltz- 
mann weights (|27[) arc independent of the representation 
of the Xi and Ei. The most illuminating representation 
of these generators is a graphical one. To find this, it is 
useful to study first another representation of the same 
algebra [H, HI ■ 

This representation arises in the Q-state Potts model 
with infinitely-strong antiferromagnetic nearest-neighbor 
interactions, and ferromagnetic next-nearest-neighbor in- 
teractions. A "spin" Si, taking the integer values Si — 
1 ... Q, is placed at each site i of the square lattice. 
The restriction on nearest-neighbor spins is simpler than 
that of the height models: adjacent spins must be dif- 
ferent. If all configurations were weighted equally, this 
would be the antiferromagnetic Q-state Potts model at 
zero temperature. However, at the critical point of inter- 
est, there are ferromagnetic next-nearest-neighbor inter- 
actions. Labeling the spins around a plaquette as Sj_i, 
Si, Sj+i and in the same fashion as the heights in 
figure [SI the restriction on nearest neighbors means that 

Pi — &Si — i Si ^SiSi-^-i fisi-i-i s' I 

for each plaquette. The transfer matrix is then given by 
with 

x l = r, (29) 

Ei ^i^Si-iSi+i 

It is easy to check that these satisfy the 50(3) BMW 
algebra with Q = (q + q^ 1 ) 2 [3^ . For Q < 4 there is a 
critical point described by a conformal field theory with 
c = 1 for Q = 3 and c = 3/2 for = 4. Using flU]) in the 
Hamiltonian ([28)) makes it clear that the next-nearest- 
neighbor interaction wants to make spins the same, so 
this critical point presumably separates an ordered phase 
from a disordered one. 

A graphical representation of /, E and X in ([29]) is 
presented in figure [T2"l The lines represent domain walls 
between different spins; any spins not separated by a do- 
main wall must be the same. The X generator has regions 
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I E X 

FIG. 12: The graphical representation of the SO(3) BMW 
algebra 

touching only at a point; these may (but need not) have 
the same spin. 

A fully packed loop model is obtained in the same fash- 
ion as the Temperley-Lieb case. The partition function 
is expanded out in powers of x and y. Each term cor- 
responds to a loop configuration after the heights are 
summed over. The generators / and E are 90-degree ro- 
tations of each other, so a rotationally-invariant model 
is obtained by taking y = 1 and x = l/(y/Q + l) -1 in 
([27]) . This yields a loop model with a typical config- 
uration illustrated in figure |U The loops here have a 
very important difference with the self- and mutually- 
avoiding loops of the O(n) loop model: because of the 
X vertex, the loops can touch. (Thus properly speaking, 
they shouldn't even be called loops, but rather graphs 
with quadrivalent vertices.) 

Since the k = 2 RSOS height models are based on a 
different representation of the same 50(3) BMW alge- 
bra, these have the same partition function as the loop 
models when y/Q = q + q~ x = 2 cos(-7r/ (p + 1)) for p inte- 
ger. Therefore the loop models at these values of p also 
have critical points with the same conformal field theory 
description, the superconformal minimal models. The 
graphical representation is valid for all Q, not just the Q 
integer where (|29|) applies, or for Q — 4cos 2 (7r/(p +1)) 
with p integer for the height models. It is likely that 
these loop models have a critical point for all Q < 4, 
although its field-theory description will only be unitary 
for p integer. 

The topological weight can be worked out using the 
BMW algebra. However, it is much simpler to find it by 
using the representation ([29]) valid for integer Q. The 
lines in the graphical representation correspond to do- 
main walls for the spins Si. This constrains spins in re- 
gions separated by domain walls to be different (if two 
regions meet only at a point, the corresponding spins are 
still permitted to be the same). Each allowed spin config- 
uration then receives the same weight, so the topological 
weight per loop configuration is given by the number of 
different spin configurations. This is simply the number 
of ways of coloring the regions with Q colors, subject to 
the constraint that adjacent regions have different colors. 
This is simply \Qi as defined in section ITT] This argu- 
ment holds true for any positive integer Q. This can be 
generalized uniquely to arbitrary Q by recognizing that 
Xq for a given dual graph satisfies the recursion relation 
©, and thus is a polynomial in Q. Thus the sum over 
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heights for any Q yields XQ- To get the full topological 
weight, note that every time an appears in the ex- 
pansion of the partition function, it comes with a weight 
x = l/(VQ + 1), as detailed after equation (|2"7]) . The 
number of X vertices Afx is a topological invariant. The 
topological weight is therefore that given in {7|) ; namely 



w(£)= X q(C)( 



1) 



This completes the mapping of the k — 2 RSOS height 
models onto the loop model described in section[H] How- 
ever, this fully packed loop model is not the only one 
in this universality class. As described in appendix [Cj 
the SO(3) BMW algebra can be written in terms of the 
Temperley-Lieb generators. In this representation, there 
are two lines on every link, with a projection operator 
Vi ensure that the two are in the spin-1 representation 
of Uq(sl(2)). Therefore, the graphical representation in 
figure [7] can be used with (|C1|) to give a fully packed loop 
model on the CuO^ lattice. The CuOi lattice is formed 
by taking a square lattice (the coppers), and adding an 
extra site on each link (the oxygens). At every copper 
site, one gets /, X, or E as above. However, instead of 
using the above graphical representation with an inter- 
section for X, one writes them instead in terms of / and 
e using (|C1|) . and uses figure [7] to represent them graph- 
ically (illustrations can be found in [32|]). The projectors 
on the links (the oxygen sites) are expanded in terms of / 
and e as well, by using Vi = 1 — ei/ (q + q~ 1 ) . This means 
that on the oxygen sites as well. This yields a somewhat 
strange but well-defined fully packed loop model on the 
CuC>2 lattice, which has a critical point described by the 
superconformal minimal models. This construction can 
readily be generalized to larger k integer as well. 



V. NEW DILUTE LOOP MODELS 

In this section I give a dilute-loop description of the 
k = 2 height models and superconformal minimal mod- 
els. The methods are similar to those used in sections 
IIII CI and IIIIDI for the k — 1 height models and the con- 
formal minimal models. As opposed to the precise results 
for the fully packed loop models, the mapping for dilute 
loops is indirect, and thus is just a conjecture, not a proof. 

The corner transfer matrix results for the fused RSOS 
models have long been known [3(| • The results of inter- 
est are along an integrable line of couplings, known as 
"Regime III" in the literature. Like for the hard-square 
and k = 1 RSOS models, along this line, the Hamiltonian 
(|28p has degenerate ground states. This is apparent from 
the computation of the local height probabilities from the 
corner transfer matrix. The formula is identical to the 
earlier case; the only distinction is that the restrictions 
on the heights are different. Namely, the local height 
probability is still given by (|19p . where q parametrizes 
the integrable line of couplings, with q — > the frozen 
limit and q — > 1 is the multicritical point. 



A. Dilute loop models for k = 2 



For k = 2, the heights are even integers hi < p, and 



hi and hj on adjacent sites obey hi 



h ; 



0,±2 and 



2 < hi + hj < 2p. It is convenient to distinguish be- 
tween heights on the two sublattices of the square lattice, 
each sublattice comprised of every other site. Denoting 
heights on one sublattice without a bar, and one with a 
bar, the adjacency diagrams for the k — 2 height models 
with p — 6 and p = 7 are given in figure 1131 Unbarred 
heights are always adjacent to barred heights, and vice 
versa. 





FIG. 13: the adjacency diagrams for the heights of the fused 
RSOS model with k = 2, p = 6 (left) and k = 2, p = 7 (right). 



The energy E in (|20j) is minimized when the heights 
on every other site are the same. This is an exact 
ground state of the model in the frozen limit. As ar- 
gued in section IIIIDI each configuration with this prop- 
erty should correspond to a ground state of the Hamilto- 
nian all along the integrable line, including at the mul- 
ticritical point, where H is given in (|28p. These ground 
states are therefore of the form G h j t = (. . . hhhh . . . ) 
and G h — (■ ■ -h,h±2,h,h±2, . . .). As before, each 
ground state corresponds to a edge on the height ad- 
jacency diagram. Excited states are defects between 
ground states. For the k — 1 case, a defect at site 
i means fti+i — hi-% — ±2. Here, a defect can have 
hi + i — hi-i = ±2 or ±4. However, for the latter de- 
fect, E in (f2T)]) is twice as large. It is therefore natural to 
assume that ±4 defects can be treated as two of the ±2 
defects. (Thus one can have a "double defect" on one pla- 
quette.) Thus fundamental defects are only those where 
hi-i — hi+i = ±2, so the adjacency diagram for ground 
states looks like those illustrated in figure [14] Another 
way of making the distinction between fundamental and 
composite defects is that a fundamental defect occurs be- 
tween G hl j t and G h2 j 2 if either hi — hi or ji = ji- 

The ground states and fundamental defects are easiest 
to list in terms of new labels. Relabel ground states of 
the form G h j l as Q\_ x , those of the form G h j^ as Q\, 
and those of the form G h+2 h as Q\. The different ground 
states Q r s then have r = 1 . . . 3, and s = 1 . . .p — 2 with 
r + s odd. The fundamental defects then separate the 
ground states G^l an d Qll \ir\ = r% ± 1 and s\ = S2 ± 1. 
Note that with this relabeling, the constraints on r and s 
are independent of each other. Moreover, the constraints 
on r and s individually are exactly the same as for the 
ground states of the k = 1 height models. This means 
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FIG. 14: The ground-state adjacency diagrams for k — 2, p — 
6 (left), and k = 2,p = 7 (right). 



that the ground-state adjacency matrix for the k 
models can be written as the tensor product 



A 2 * = A 1A <g> A hp - 



(30) 



where we denote the adjacency matrices for general k and 
p by A k,p . The adjacency matrix for the heights does not 
obey such a simple decomposition: this is only a property 
of the ground states and fundamental defects. 

As with the k = 1 case, the world lines of the fun- 
damental defects are interpreted in the two-dimensional 
lattice model as domain walls between free-energy min- 
ima. These domain walls form the loops. Domain walls in 
these RSOS models occur when heights across a square 
plaquette are different. These domain walls are drawn 
on the diagonals of the original lattice, as seen in the 
(non-intersecting) loops of figure [9] (In the fully packed 
loop models, the links are drawn on the links of the dual 
lattice.) Each link of the original lattice corresponds to 
a ground state, so the domain walls separate links with 
different ground states. 

To complete the argument, the topological weight w(C) 
of each loop configuration must be found. The simplest 
possibility, realized for k = 1, is that is that these do- 
main walls never cross, so that the loops self-avoid and 
mutually avoid. Summing over the heights in the model 
then gives the topological weight of the 0(n) loop model: 
w = , where Af is the number of loops. For several 
reasons, such a w cannot be correct for k > 1. 

First of all, domain walls do cross for k > 1. This 
happens when all four ground states around a link are 
different. In the k = 2 model, the domain walls cross 
when the heights around a plaquette are of the forms 



h + 2 



h 




h + 2 




h + 2 



h + 2 



(31) 



The four different ground states on these two plaquettes 

are G h,h > G h+2,h > G h+2,h+2 &nd G h,h+2 ■ These have 
a non-zero Boltzmann weight coming from Xi, as given 

in appendix [C] Configurations like these cannot occur 

in the k = 1 model, where ij requires that the top and 

bottom heights be the same, while requires that the 

left and right be the same. 



Even if one were to assume that such crossings are 
irrelevant, the value of n obtained is still inconsistent. 
Namely, if crossings are ignored, one obtains the dilute 
0(n) loop models as with k — 1. The arguments of sec- 
tion IIII Dl give the weight n per loop to be the largest 
eigenvalue of A. This is easy to find from (f30|) : it is 
simply the product of the largest eigenvalues of A 1 ' 4 and 
A 1 *- 1 , which are 2cos(tt/4) = y/2 and 2cos(?r/(p- 1)) 
respectively. This would yield n = 2\p2 cos(ir/(p — 1)). 
Although this is the precise value of n implied by the SLE 
results [9(, the loop model cannot be the 0(n) model. 
The weight per loop has n > 2 (e.g. n — 2\[2 for the 
5?7(2)2 WZW model) in general, and the 0(n) model in 
this regime does not have a critical point. Since the un- 
derlying lattice model is critical, a non-critical loop model 
obviously cannot be an equivalent description. The fact 
that domain walls can cross or touch must affect the con- 
tinuum behavior of the loop model. 

The fact that the weight of non-crossing loops is the 
same as the SLE results is quite encouraging. The task is 
then to use the Boltzmann weights of the height model to 
understand how to treat the crossings. The fully packed 
loops do intersect, and affect the topological weight via 
the chromatic polynomial. The dilute loops behave dif- 
ferently. The decomposition (f3"0"l) suggests that for k > 1 
the domain walls should not be treated as a single lines, 
but instead split into two lines, each of which behaves as 
a k = 1 model. 

The key to defining a topological weight is the rela- 
beling of the ground states as as described above. In 
this relabeling, the r and s labels behave independently, 
as far as the adjacency restrictions on the ground states 
are concerned. Since the domain walls are between the 
ground states, this means the domain walls in the r and 
s labels can be treated independently. Thus the domain 
walls can be split into a "solid" wall and a "dashed" wall, 
as illustrated in figure [21 Away from where the domain 
walls meet, this splitting is unimportant, and the solid 
and dashed lines are attached. However, at a domain- 
wall crossing there are now four ways of resolving the 
lines so that lines of the same type do not cross. These 
are illustrated in figure 1151 The dashed line is a domain 



J 1 ... 



FIG. 15: The four ways of resolving a domain- wall crossing 
with doubled lines 

wall in the 1,4 theory (the r index in Q r s ) : while the solid 
line corresponds to a domain wall in the l,p — 1 model 
(the s index in Q r s ). 

Note that dashed lines never cross or touch each other, 
and likewise for solid, but that the dashed and solid can 
cross one another. Thus a loop configuration consists of 
dashed and solid loops which are sewn together, except 
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when two loops meet on one plaquette. At each domain- 
wall crossing, one of the four possibilities in figure 1151 
must occur. A typical loop configuration is illustrated 
in figure O Since the restrictions on the r and s labels 
are completely independent, it is consistent to treat the 
dashed and solid domain walls independently. 

All the configurations in the height model can therefore 
be realized in terms of a pair of k = 1 domain walls. 
The last two domain- wall configurations in figure [TBI are 
necessary to realize configurations like (|31[) that occur 
for k > 1. This follows from rewriting the ground states 
around the plaquette in terms of Q r s . For example, G h - t — 

Q\_ x and G h+2 ^p2 = Qh+i- Therefore they have the 
same r value and are not separated by any dashed lines. 
Likewise G h+2 - t = Q\ and G h = G\ so these have the 
same s value and are not separated by any solid lines. 

This picture is in remarkable accord with the field- 
theory results for the (k = 2) superconformal minimal 
models perturbed along the integrable line of interest. 
The particles in the field theory have the same labels as 
the defects, so the particle world-lines can again be iden- 
tified with the domain walls. Moreover, the scattering 
matrix S 2 ' p is a product of k = 1 scattering matrices 
[H]], just like the ground-state adjacency matrices, ft is 

S 2 'P = S 1 ' 4 ® S 1 *- 1 . (32) 

Plugging the k = 1 scattering matrices from (|23|) into 
(f3"2"| shows that S 2p is a sum of four terms, involving 
/(g)/, e ® /, / ® e, and e <g e. Recalling the graphical 
representation of / and e in figure [JJ gives precisely the 
four possibilities in figure [15] for the particle world-lines. 

It is now clear how to build a topological weight for 
these dilute loops. The dashed and solid lines each form 
closed loops. The sum over heights can be done indepen- 
dently for the r and s labels for each ground state. The 
weight per dashed loop is then simply the largest eigen- 
value of the A 1 ' 4 adjacency matrix, which is \/2. The 
weight per solid loop is the largest eigenvalue of A 1 ' 1 '" 1 , 
which is 2cos(7r/(p — 1)). When the numbers of dashed 
and solid loops are M. and M respectively, the topological 
weight is 

w{C) = V2 M (2 cos (— rr)) b ° ■ ( 33 ) 

The topological invariant C is the number of times the 
(original) domain walls cross, i.e. the number of plaque- 
ttes which have one of the configurations in figure [15] at 
their center. The piece involving the parameter b is anal- 
ogous to the last piece of the fully packed weight ([7]). The 
arguments here are not refined enough to determine the 
value of 6, or whether changing it results in a relevant or 
irrelevant perturbation of the critical point. 

B. Arbitrary k 

The ground states for general k are quite analogous to 
those for k = 1 and k = 2. The expression (|19l20p for 



the local height probabilities still holds, so the minimum 
of E occurs when the heights on every other site are the 
same. Thus each ground state can be labeled by two 
heights hi and hj, which must satisfy the constraints for 
nearest- neighbor heights in ([26]) . 

For general k and p integer, the adjacency matrix de- 
composes into tensor products of k = 1 adjacency matrix. 
This follows by relabeling the ground states Gi n ,h 2 into 
the form Q r s via 

hi - h 2 + k + 2 

T = 

2 

hi + h- 2 — k 



The constraints in (|26p in terms of r and s are therefore 
quite simply r = 1, 2 . . .k + 1 and s — 1, 2, . . .p — k. 

A fundamental defect needs to have energy E = 2 in 
([20]) . This only occurs if the two ground states have 
one of their height labels in common, and one different. 
Namely, a fundamental defect separates the ground state 
G hl ,h 2 lrom the ground states Gh 1 ,h 2 ±2 and G hl ±2,h 2 - In 
terms of r and s labels, this means that a fundamen- 
tal defect separates Q r s from the states Gl±l and G^±l- 
Thus the restrictions on r and s are independent of each 
other, and the adjacency matrix for the ground states for 
general k and p is simply 

A k, P = A i,k+2 9 A i,P-k+i ( 34 ) 

This is in harmony with the scattering matrix for the 
quasiparticles in the field theory describing the integrable 
line, which is [35] 

S k * = S 1 > k+i ®S 1 *- k+1 . (35) 

This also is in harmony with the Coulomb-gas approach 
to these conformal field theories, which gives their par- 
tition functions as sums of products of minimal-model 
ones [Hj]. 

The configurations in the dilute loop model for gen- 
eral k are therefore the same as for k = 2 model. The 
defect world-lines become "double" domain walls in the 
two-dimensional picture. Crossings are resolved in the 
four ways in figure (|15p . so the typical configuration dis- 
played in figure [2] still is applicable here. The weight for 
the dashed loops for general k is now the largest eigen- 
value of A 1 '^ 2 , which is 2cos(7r/(fc + 2)). This yields the 
topological weights given in (|3l4p . 

VI. FURTHER DIRECTIONS 

In this paper I have shown how to turn lattice models 
into loop models. Dense loop models arc found directly 
by exploiting the algebraic structure of the transfer ma- 
trix, while dilute loops are the domain walls found more 
indirectly using the corner transfer matrix and the scat- 
tering matrix. Since these lattice models have critical 



17 



points described by known coset conformal field theories, 
this means that the loop models have these same critical 
points as well. 

It would be interesting to show directly that the dilute 
loop models have these critical points. One way of doing 
so would be to study the "dilute BMW models" [37j [. 
These lattice models generalize the dilute Temperley- 
Lieb models [2(| Hl| discussed above to allow a more 
general algebra. The domain walls for the dilute 5*0(4) 
BMW model are exactly the same doubled domain walls 
as in the dilute loop model discussed here (because the 
algebra 50(4) = SU(2) x SU(2)). The corner transfer 
matrix computation for these models has not yet been 
done, so it is not yet known if these models have the 
critical points described in this paper. 

It would also of course be interesting to generalize these 
results to other loop models and conformal field theories. 
Another loop model with a critical point described by the 
conformal minimal models was discussed at length in [l2| . 
This dilute loop model can be derived from the Potts lat- 
tice model from both the domain-wall and the algebraic 
approach (using projection operators instead of /, X and 
E) , but the heuristic approach used there is more sugges- 
tive. This is the scattering matrix approach used above, 
where the world lines of the particles in the field the- 
ory are identified with the loops. As mentioned briefly 
in the scattering matrix approach suggests that the 
dilute version of the fully packed k — 2 model discussed 
above will have critical points described by parafermion 
conformal field theories. This could be tested directly by 
studying the 50(3) dilute BMW model of 

Many integrable lattice models and scattering matrices 
are known, and it likely that the methods described here 
could be applied to them. Finding more examples would 
be of great interest in both the areas mentioned in the 
introduction: SLE and non-abelian statistics. Much is 
known about the relations between conformal field the- 
ory, non-abelian statistics, topology and geometry in the 
simplest cases. Hopefully the generalizations described 
here will be of use in deepening this knowledge. 
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APPENDIX A: A BRIEF REVIEW OF THE 
CONFORMAL FIELD THEORIES 

Field theories describing rotationally-symmetric two- 
dimensional classical critical points are invariant un- 
der conformal transformations of two dimensional space. 
Conformal symmetry in two dimensions is very powerful, 
because it has an infinite number of generators: represen- 
tations of the corresponding symmetry algebra (called 



the Virasoro algebra) are infinite-dimensional. Mini- 
mal models are conformal field theories which have a fi- 
nite number of highest- weight states under the conformal 
symmetry These are among the models for which 
SLE is applicable. The minimal models have been clas- 
sified, and all the critical exponents are known 0, l38| . 
Many (and in principle, all) of the correlation functions 
can be explicitly computed. A convenient way of labeling 
conformal field theories is in terms of a number c called 
the central charge. Unitary minimal models have central 
charge 



G 



p(p+l) 



(Al) 



where p is an integer obeying p > 3. These are the only 
unitary conformal field theories with c < 1 [38]. Any 
theory with c > 1 has an infinite number of highest- 
weight states (39j . 

To study more general conformal field theories, it is 
useful to extend the Virasoro algebra by some other 
generators. One well-studied way of doing so is to ex- 
tend the symmetry by a simple Lie algebra G. These 
conformal field theories are called Wess-Zumino-Witten 
( WZW) models [40] , and the extended symmetry algebra 
is called a Kac-Moody algebra. Like with the minimal 
models, all critical exponents are known, and correlators 
can be computed explicitly 41]. Each unitary WZW 
model Gk is labeled by the algebra G and an integer 
k > 1, called the level. 

Minimal models can be constructed from WZW mod- 
els by using the coset construction [42^ . A G/H coset 
conformal field theory is defined with energy-momentum 
tensor T G m = Tg~Th, where Tq and Tjj are the energy- 
momentum tensors of the G and H WZW models. This 
allows many of the properties (e.g. the critical exponents) 
of the coset models to be computed by using results from 
the WZW models. The coset models of interest here are 



SU(2) k x SU(2) P ^ 1 
5C/(2) P _! 

with p > k + 1. They have central charge 



3k 
k + 2 



1 - 



2(fc + 2) 



(p + l)(p-*E + l) 



(A2) 



(A3) 



When k = 1, these are the minimal models. The mod- 
els with k — 2 are usually known as the Af — (1,1) 
superconformal minimal models, because they have an 
extended symmetry algebra including supersymmetry. 
Note that the simplest superconformal minimal model 
(with k — 2,p = 4) is identical to the second confor- 
mal minimal model (with fc = 1, p = 4). As p — > 00 
for fixed k, one obtains simply the SU(2)k WZW model. 
The Coulomb-gas formulation of these coset models was 
developed in [36j . 

In the loop models discussed below, p becomes a con- 
tinuous parameter. The conformal field theories describ- 
ing loop models for p non-integer generally are expected 
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to be non-unitary and not rational. Exact computations 
are still possible by using the Coulomb-gas technique 
within conformal field theory, and (for c < 1) SLE. 

APPENDIX B: LOOPS ON THE SPHERE 

All the results of this paper apply when two- 
dimensional space is topologically a sphere. In this ap- 
pendix I show how to implement these boundary con- 
ditions with the transfer matrix written in terms of 
Temperley-Lieb generators, and show that the partition 
functions in this situation are independent of representa- 
tion of this algebra. 

In a loop model, space is topologically a sphere when 
no loops end on the boundary. To implement these 
boundary conditions in the transfer-matrix formulation 
of the height model, define the operator 

B = (q + q~ x ) N eies ... eajv-i, 

which acts on a zig-zag row of heights as illustrated in 
figure [H acting on a set of heights h®, h\, . . . h 2 N- B is 
normalized so that B 2 = B. Then the partition function 
for M + 1 zig-zag rows is 

Zheight = (ST S){ a }{ a } 

where {a}{a} means the matrix element which has all 
the heights in the first and last rows fixed to be a, a + 
1, a, a + 1, . . . . The fact that T involves only Jo and I 2 n 
for the left- and right-most columns means that all the 
heights along the sides are fixed to the same values. 

To relate this definition of Z with that coming from 
the loops, expand T in terms of the Temperley-Lieb gen- 
erators, as done for the fully packed loop models. Each 
of the resulting terms corresponds to a loop configuration 
£, and consists of many sandwiched in between the Bs. 
Consider one of these terms, which gets weight A(C) in 
the loop model, and A^ a j^ a y in the height model. Then 
using the algebra (fT5]) , I show here that 

A(C)B {a}{a} = A {a}{a} . (Bl) 

When computing A(C) the in B are drawn with only 
one loop: the piece that goes "off the edge" is ignored. 
The relation (|B1[) shows that up to an overall constant, 
the partition functions are the same: 

z loop = S MM z height ( B2 ) 

Let j be an integer. Because (e;) 2 = (3 + q~ l )e.i and 
all commute with each other, 

e 2j -iB = Be 2j -i = {q + q- l )B. 

Thus any time an e 2 j— 1 is next to B, it can be removed, 
leaving an overall factor q + q^ 1 . Now consider an e 2 
somewhere in the middle. All one needs to do to get 
rid of it is to get the eis from the Bs next to the e 2 , and 



then use eie 2 ei = e±. Since e\ commutes with everything 
other than e 2 , the only obstruction to doing this is the 
product is of the form 

eie 2 . . . e 3 . . . e 2 ei. 

where e\ does not appear in the dots. Since there are no 
e]S in the dots, e 2 can be commuted through the dots 
and e 2 e 3 e 2 = e 2 used, unless there is an e 4 within the 
dots, i.e. the product is of the form 

eie 2 e 3 . . . e 4 . . . e 3 e 2 ei. 

The offending e 4 can then be eliminated using the e 3 s 
unless there is an e§ which interferes. Thus either the 
offending ones can be eliminated, or one ends up with 

eie 2 e 3 e 4 . . . e 2N - 2 e 2N -ie 2N - 2 . . . e 3 e 2 ei. 

But then this can be reduced by using 
e 2 Ar- 2 e 2 Ar_ie 2 Ar_ 2 = e 2A r_ 2 and so on. Proceeding 
like this one can eliminate all of them, leaving 

A=(q + q- 1 fl3 (B3) 

where Af is an integer independent of representation. In 
fact, Af is the number of closed loops in the graphical rep- 
resentation of C. This proves (|B1|) and (|B2p for any rep- 
resentation of the ej, as long as it satisfies the Temperley- 
Lieb algebra. 

For the loop models based on the 50(3) BMW algebra, 
the e in B are replaced with E. The arguments then go 
through as above. 

APPENDIX C: E, X AND BMW 

This appendix collects formulas connected to the 
SO(3) BMW algebra. Much of this appendix is discussed 
in depth in [32j , and what is not can be obtained by using 
these results along with the exp licit Boltzmann weights 
for the k = 2 RSOS models [if. 

The generators Xi and Ei can be defined in terms of 
the Temperley-Lieb generators by using the projectors 
Vi = I — ei/(q + cT 1 ) onto the spin-1 representation of 
U q (sl(2)), giving |H 

Ei = V 2 i+iV 2 i+3 e 2 i+ 2 e 2 i + i e 2 i + 3 e 2 i +2 V 2 i+iV 2 i+3, 

Xi = {q + q- 1 )V2i+ir2i+3e 2i+2 V2i+iV2i+3- (CI) 

Representations of the SO(3) BMW algebra therefore 
can be found from representation of the Temperley-Lieb 
algebra simply by plugging into (|C1[) . In particular, a 
graphical representation in terms of non-crossing loops 
can be found from the pictures in figure [7] 32]. 

The S'0(3) BMW algebra is usually written in terms of 
Ei and -B^, where Bi is a braiding operator (i.e. satisfies 
BiBi + \Bi = Bi + iBiBi + i). Here it is convenient to write 
them in terms of I,X, and E; the expression for B in 
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terms of these can be found in [32j . The Ei satisfy the 
Temperley-Lieb algebra amongst themselves: 



Ef = {Q-l)E u 
EiEi±\Ei = Ei. 

Relations involving the X are: 

(Xi) 2 = (Q-2)X l + E l 



(C2) 
(C3) 



and 

XiEi + iX.- L 
EiXi + iEi 
X{Ei + iEi 
XiXi + iEi 
XiXi + iXi — Xi + iXiXi + i 



Xi+iEiXi + i, 

(Q-i)E h 

Xi+iEi, 

(Q — 2)X i+ iEi + Ei, 
Xi+\Ei + EiXi + \ + Xi 
—Ei — XiEi + \ — Ei + iXi 
—Ei+i + Xi + i. 



All relations also hold with i and i + 1 interchanged, and 
with the order of products in each term reversed. All 
generators labeled by sites i and j commute when \i— j\ > 
1 . It is instructive to use the graphical representation in 
figure [T2l to verify these relations. 

Explicit expressions for E and X in the k = 2 height 
representation discussed in this paper can be found either 
by substituting (fTT|) into (jClj) , or from the explicit Boltz- 
mann weights. Since the Ei satisfy the Temperley-Lieb 
algebra, they end up satisfying the identical relation as 
the ei , the only difference being that the heights here are 
always even, with the adjacency rules in ([26]) for k = 2. 
Define [a] = sin(a7r/(p+ 1)). The X, in the height repre- 



sentation are 

h + 2 






h + 2 





[h+1] 



[h+3] 
[h+1] 



[h+m\ 
[h+i][i] 



[h-l] / [h+2] 

[h+i]V [h] 



ji]_ (\h-if | [h+i] 2 



[2}[h] V [h+1] 1 [h-l] 



The remaining non-zero elements of Xi are found by ei- 
ther using the symmetry h <-> p + 1 — h, or by flipping 
them horizontally or vertically. All others are zero. The 
fact that the first of these is non- vanishing has a number 
of interesting consequences discussed in section IVl 
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